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ABSTRACT 

Previous work on the sonic -point beat-frequency (SPBF) model of the kilohertz quasi-periodic oscillations 
(QPOs) observed in the X-ray flux from neutron stars in low-mass binary systems has shown that it naturally 
explains many properties of these QPOs. These include the existence of just two principal QPOs in a given 
source, the commensurability of the frequency separation Aiy of the two kilohertz QPOs and the spin frequency 

inferred from burst oscillations, and the high frequencies, coherence, and amplitudes of these QPOs. Here we 
show that the SPBF model predicts that Af is less than but close to fs, consistent with the observed differences 
between Af and Vs- It also explains naturally the decrease in Aiy with increasing QPO frequency seen in some 
sources and the plateau in the QPO frequency-X-ray flux observed in 4U 1820—30. The model fits well the QPO 
frequency behavior observed in ScoX-1, 4U 1608-52, 4U 1728-34, and 4U 1820-30 (xVdof = 0.4-2.1), 
giving masses ranging from 1.59 to 2.Q Mq and spin rates ranging from 279 to 364 Hz. In the SPBF model, the 
kilohertz QPOs are effects of strong-field gravity. Thus, if the model is validated, the kilohertz QPOs can be used 
not only to determine the properties of neutron stars but also to explore quantitatively general relativistic effects 
in the strong-field regime. 

Subject headings: relativity: general — stars: neutron — stars: oscillations — stars: rotation — stars: X-ray 



L INTRODUCTION 

Kilohertz quasi-periodic brightness oscillations (QPOs) have 
been discovered in the accretion-powered emission of more 
than twenty neutron star low-mass X-ray binaries using the 
Rossi X-ray Timing Explorer (RXTE) (see van der Klis 2000 
for a review). These kilohertz QPOs have frequencies ranging 
from ~500 Hz up to '-1200 Hz, ampHtudes --1%-15% rms in 
the 2-60 keV band of the RXTE Proportional Counter Array, 
and quality factors Q = v/FWBM as lai-ge as 100-200. 
They often appear as two simultaneous peaks in power spec- 
tra of the X-ray brightness, at frequencies ui and 1^2 (> i^i) that 
move together and can shift upward and downward in frequency 
by up to a factor of two within a few hundred seconds, appar- 
ently as a consequence of changes in the accretion rate. The 
frequencies of the kilohertz QPOs lie within the range expected 
for orbital motion near neutron stars and are almost certainly 
strong-field general relativistic phenomenon (see Lamb, Miller, 
& Psaltis 1998; Miller 2000; van der KHs 2000). 

In six of the kilohertz QPO sources, strong (amplitudes up 
to 50% rms) brightness oscillations with frequencies of hun- 
dreds of Hertz have been observed during thermonuclear X-ray 
bursts. There is very strong evidence that these oscillations are 
produced by rotation of brighter regions of the stellar surface 
with the spin of the star and that the most prominent oscillation 
frequency is the stellar spin frequency or its first overtone 
(see, e.g., Strohmayer & Markwardt 1999; Miller 1999). The 
coherent oscillations observed in burst tails, after hot matter 
has spread around the star, and the antipodal brighter regions 
observed in 4U 1636—53 are strong evidence that these stars 
have dynamically important surface magnetic fields (Miller 
1999), as expected on evolutionary grounds (Miller, Lamb, & 



Psaltis 1998, hereafter MLP98; Lorimer 2000). In the four 
sources where burst oscillations and two simultaneous kilohertz 
QPOs have both been detected with high confidence, the differ- 
ence Af = V2 — vi between the frequencies of the kilohertz 
QPOs agrees with the fundamental burst oscillation frequency 
to within 0.7% to 15% (see van der Klis 2000). 

The properties of the kilohertz QPOs strongly indicate that 
the upper kilohertz QPO is an orbital frequency and that the 
lower is generated by the beat of this frequency with the star's 
spin (Lamb et al. 1998; MLP98; Miller, Lamb, & Cook 1998): 
( 1 ) The frequencies of the kilohertz QPOs are in the range ex- 
pected for orbital motion near neutron stars, consistent with one 
being an orbital frequency. (2) The relatively small variation of 
Av in a given source and the approximate commensurability 
of Av with the stellar spin frequency inferred from burst os- 
cillations indicates that the spin of the star is generating the 
frequency difference. (3) The observation of at most two strong 
kilohertz QPOs in any given source (van der Klis 2000; Mendez 
& van der Klis 2000), rather than three, indicates that only a sin- 
gle sideband is being generated, rather than the two sidebands 
generated by other mechanisms. (4) In order to beat with the 
stellar spin, the other motion must be rotation about the center 
of the neutron star, indicating that the other motion is orbital 
motion. (5) Accretion of angular momentum ensures that the 
star is spinning in the prograde direction, implying that the fre- 
quency V2 of the upper kilohertz QPO is close to the orbital 
frequency t^orb and that the frequency vi of the lower kilohertz 
QPO is close to the beat frequency vb = Voih — i^s- 

This and other evidence motivated the development of the 
sonic-point beat-frequency (SPBF) model (MLP98). In this 
model ^2 is close to the orbital frequency I'orb at the sonic ra- 
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dius, where the flow in the disk changes from nearly circular 
to rapidly inspiraling, and i^i is comparable to i^b- The tran- 
sition to hypersonic radial inflow is an effect of strong-field 
gravity. In the SPBF model, the low-frequency QPOs seen 
in the Z and atoll sources are produced by the magnetospheric 
beat- frequency mechanism (Alpar & Shaham 1985; Lamb et al. 
1985). The model explains naturally many features of the kilo- 
hertz QPOs in addition to their frequencies and frequency sep- 
aration and the existence of only two principal QPOs in a given 
source (MLP98). In the SPBF model, the kilohertz QPOs are 
a strong-field general relativistic effect and are therefore sensi- 
tive probes of the properties of the spacetime near the neutron 
star, including whether there is an innermost stable circular or- 
bit (ISCO) and the gravitomagnetic torque created by the spin 
of the neutron star. 

Careful analyses have shown that in some kilohertz QPO 
sources Ai' decreases systematically by 30-100 Hz, depend- 
ing on the source, as 1^2 increases by a much larger amount 
(Sco X-1: van der KHs et al. 1997; 4U 1608-52: Mendez et al. 
1998; 4U 1735-44: Ford et al. 1998; 4U 1728-34: Mendez & 
van der Klis 1999; see also Psaltis et al. 1998) and that the sep- 
aration frequency in 4U 1636—536 is slightly but significantly 
smaller than the spin frequency inferred from its burst oscilla- 
tions (Mendez, van der Klis, & van Paradijs 1998). 

Here we show that in the SPBF model, the inward radial ve- 
locity of the accretion flow increases and decreases 1^2, mak- 
ing Aiy smaller than h's- An increase in vi of <^ 5% is suffi- 
cient to account for the largest observed differences between 
Al/ and i^b and decreases of Ai' with increasing 1/2- The in- 
ward radial velocity of the accretion flow was included in the 
gas dynamical and radiation transport calculations reported by 
MLP98, but the effect of this motion on the frequencies of the 
kilohertz QPOs was neglected. In § 2 we discuss the generation 
of the kilohertz QPOs in the SPBF model, derive general ex- 
pressions for their frequencies, taking into account the inward 
motion of the accretion flow, summarize the results of numer- 
ical simulations, and introduce the model of the gas dynamics 
in the SPBF model that we use here to calculate the frequencies 
of the kilohertz QPOs. In § 3 we show that despite its simplic- 
ity, this model fits accurately the frequency behavior observed 
in Sco X-1, 4U 1608-52, 4U 1728-34, and 4U 1820-30 and 
gives reasonable masses and spin rates for these neutron stars. 

2. KILOHERTZ QPOS IN THE SONIC-POINT MODEL 

Generation of kilohertz QPOs. — In the SPBF model, the 
kilohertz QPOs are produced by clumping of gas in the accre- 
tion flow near the star, probably caused by interaction of the ac- 
creting gas with the weak stellar magnetic field (see Ghosh & 
Lamb 1979; Aly & Kuipers 1990). Only clumps near the sonic 
radius i?sp, where the inward radial velocity of gas leaving the 
clumps becomes supersonic, contribute to the oscillations (see 
MLP98). At low accretion rates, radiation drag removes angu- 
lar momentum from the gas in the clumps, causing it to spiral 
inward to the stellar surface. At high enough accretion rates, 
the optical depth from the stellar surface is large and radiation 
drag is unimportant. If, however, the star is smaller than the 
ISCO, gas will fall inward supersonically from clumps near the 
ISCO, where general relativity causes gas to spiral inward even 
if it does not lose any angular momentum. In either case, the 
supersonic inflow is a strong-field general relativistic effect. 

Where gas falling inward from clumps near R^-p impacts the 
stellar surface, it produces bright footprints that move around 
the star with a frequency nearly equal to the orbital frequency 



forb of the clumps at i?sp, causing the X-ray flux seen by dis- 
tant observers to oscillate as the footprint moves into and out 
of view. (As discussed by MLP98, i^orb is not exactly equal to 
the Keplerian frequency vk, because the radiation force has a 
radial component and the clumps are drifting slowly inward.) 
This oscillation is the upper kilohertz QPO; its frequency V2 
is determined by the rate at which the azimuthal phase of the 
footprint advances. 

The weak magnetic field of the star funnels extra gas toward 
the magnetic poles, producing a broad beam of extra radia- 
tion that rotates with the star. Clumps orbiting near Rgp move 
through this beam with frequency i^b = ^'oib — v^. When a 
clump is illuminated by the beam, gas in it loses angular mo- 
mentum to the radiation at a faster rate, so more gas falls in- 
ward, causing the luminosity of the footprint to oscillate with 
a frequency close to vb- This oscillation is the lower kilohertz 
QPO; its frequency vi is determined by the frequency at which 
the massfiux onto a footprint oscillates. 

The flux in both oscillations is generated primarily at the stel- 
lar surface, where most of the gravitational energy is released, 
and hence both oscillations can have relatively high amplitudes. 

Kilohertz QPO frequencies and frequency separation. — The 
kilohertz QPO frequencies predicted by the SPBF model can 
be calculated as follows (we use Boyer-Lindquist coordinates 
for easy reference to infinity). The light travel time from the 
stellar surface to the sonic radius is very small and hence the 
time between two successive maxima of the mass inflow rate 
from a given clump is very close to the beat period AIb = 
l/(t'orb ~ t^s)- Let ti and ^2=^1+ AtB be the times of suc- 
cessive maxima, let ri and r2 and 2TTipi and 2TTip2 be the radial 
and azimuthal coordinates of the clump at ti and t2, let Ati and 
At2 be the times required for the gas stripped from the clump 
at ti and t2 to reach the surface of the star, and let 27rA0]^ 
and 27rA(/)2 be the changes in the azimuthal phase of the gas 
released at and ^2 as it falls from the clump to the stellar 
surface. Then the frequencies of the lower and upper kilohertz 
QPOs are 



t2-ti + {At2 - Ah) 

and 

_ 02-</'l + (A'/'2-A0i) 

t2-h + {At2 - Ah) 

Suppose first that the inward motion of the clumps near Rgp 
can be neglected. This is the approximation used in MLP98 
to compute the frequencies of the kilohertz QPOs. Then the 
time required for gas to spiral inward from a given clump to 
the stellar surface and the change in the azimuthal phase of the 
gas during its inspiral do not depend on when the gas separated 
from the clump, i.e., At2 = Ati and A(j)2 = A(/)]^. Because 
A<2 = Ail, the time between two successive maxima of the 
mass flow rate at the stellar surface is the same as at the clump, 
where it is AtB, and hence the frequency I'l of the luminosity 
oscillation is I /AtB = t'orb — J^s = i^b (see eq. [|l]]). Be- 
cause the rate at which the azimuthal phase of the footprint ad- 
vances is the same as the rate at which the azimuthal phase of 
the clump advances, the rotation frequency 1/2 of the footprint 
is {4>2 — 0i)/(^2 — ^i) = t'oib (see eq. [^]) and the difference 
between 1^2 and is therefore Vg- 

Suppose now that the inward motion of the clumps near Rgp 
is too large to be neglected. The radial distance of a given clump 
decreases with time and therefore the time required for gas to 
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spiral inward from the clump to the stellar surface steadily de- 
creases. Hence At2 < A<i. Equation (|l]) shows that to first 
order in Vd, the frequency of the lower kilohertz QPO is 
+ VcidrAt) > vb; vi is greater than vb because the in- 
ward motion of the clump reduces the spatial separation be- 
tween successive maxima of the mass flow rate, analogous to 
the Doppler shift of a sound wave. Equation (H) shows that to 
first order in v^x, the frequency V2 of the upper kilohertz QPO 
is Vorh{^ + VcidrAt) — Vc\drA<j> — i^B ^ ^ Wci^r i^orb- Numerical 
simulations of the gas flow and radiation transport in full gen- 
eral relativity similar to those presented in MLP98 show that the 
term proportional to 9^ A0 is generally the dominant correction 
term. This term describes the effect on 1^2 of the fact that the gas 
falling inward from the clump winds a progressively smaller 
distance around the star as the clump drifts inward, because the 
gas falls through a smaller radial distance before colliding with 
the stellar surface. Hence 1^2 is generally less than I'oib- The 
fractional decrease in 1^2 is typically ~ 50% of the fractional 
increase in I'l. The increase in I'l and the decrease in 1^2 reduce 
the frequency separation Aiy, which is therefore less than i^s- 

Inward motion of clumps and gas. — The previous analysis 
shows that Ai/ is less than i/g. A model of the inward drift of 
clumps near R^p and of the infall of gas from the clumps is re- 
quired in order to compute the detailed behavior of i^i and 1^2 . 
Here we consider a simple model that captures the basic physics 
involved. 

Let Vc\ be the inward radial velocity of clumps near R^p, 
let Vg be the characteristic initial inward radial velocity of the 
gas stripped from a clump, and assume that Vc\ ^ Vg. Then 
vi ^ VB/i^ — Vci/vg) (seeeq. [|l|]), because the motion near the 
clump contributes most to At. Based on the results of the nu- 
merical simulations cited above, we assume that the fractional 
decrease in 1^2 is equal to one-half the fractional increase in i'l. 
Then ~ t'oib(l — ivc\/vg). Our simulations indicate that the 
initial inward radial velocity of stripped gas is fairly insensitive 
to the distance of the clump from the star and hence for simplic- 
ity we assume that Vg is independent of R^p. We also assume 
for simplicity that Vc\ and Vg are approximately constant during 
the lifetime of a clump. 

The inward radial velocity of a clump with angular momen- 
tum J caused by removal of its angular momentum by a brak- 
ing torque N is Vc\ = N/{dJ/dr). For a clump of mass 
m in a nearly circular orbit in the Schwarzschild spacetime, 
dJ/dr = {m/2){r - 6M){r - 3Af)-3/2 ^j^^ j^gj^^g 



Wcl(^sp) = 



(2/m) (R,p - 3M) 



3/2 



{Rsp - 6M) 



-N 



(3) 



in units with G = c — 1. This expression for Vc\ diverges at 
Rsp — 6M because the inertia of the clump has been neglected. 
In reality, Vci(i?sp) increases as Rsp approaches the radius of 
the ISCO but remains finite (see Fig. |l|). 

There are several effects that may produce a braking torque 
N on orbiting clumps. Here we consider one likely possi- 
biUty. Rsp lies inside the star's magnetosphere (MLP98), 
where the pressure of the stellar magnetic field compresses 
gas clumps, isolating them from one another We there- 
fore assume that clumps near Rsp are in pressure equilibrium 
with the stellar magnetic field, that they interact only weakly 
with each other, and that N is dominated by the interaction 
of the clump with the stellar magnetic field (see Ghosh& 
Lamb 1979; Aly & Kuipers 1990; Ghosh & Lamb 1991). 




6 6.5 7 7.5 6 

Sonic Radius (M) 

Fig. 1 . — Inward radial velocity v^i of clumps orbiting near the sonic radius 
Rsp as a function of Rsp, computed numerically (solid lines) and using equa- 
tion (dotted lines), for (top to bottom) cjy = 10"^, 10"^, and 10^'^ in 
units of the star's gravitational mass M. 

Then N = ^ira^RspBsp, where 7 is a dimensionless factor 
where 7 is a dimensionless factor ^ 1 that depends on the 
distortion of the stellar magnetic field by the clump and the 
clump's shape, a is the mean radius of the clump, and Bsp is 
the magnetic field at Rsp. Radiation from the star is likely to 
maintain clumps at the local Compton temperature, which is 
almost independent of the distance to the star; we therefore as- 
sume that the clump temperature is approximately independent 
of Rsp. For simplicity, we also assume that m is approximately 
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independent of i?sp. Then {S/A'K)2mkBT/r 
H^R^p for a dipolar field with moment /z, where nip is the mass 
of a proton, and hence a oc R^p and N oc R~p. With these as- 
sumptions, the evolution equation for the specific angular mo- 
mentum Ucf, of a clump is du^/dr — cpf/r, where the torque 
coefficient c n (which like and the proper time of the clump 
has the dimensions of a mass) depends on geometrical factors, 
fundamental constants, and 7. The evolution equation for the 
radial position of a clump moving in the Schwarzschild space- 
time in the absence of significant nongravitational radial forces 



dr 



(4) 



We computed fci(i?sp) by integrating numerically the clump 
equations of motion, starting with the values of f and f given 
by equation (^ at large Rsp. 

Figure [l| shows how Vd varies with Rsp for c^v = 10^^, 
10^^, and lO^'^ in units of the star's gravitational mass M. 
These values correspond to clump densities ~ 1-10 times 
the mean density in the disk, for Bsp ^ 10^-10^ G. The in- 
ward velocity of clumps increases as Rsp approaches the ISCO 
at 6M, but remains finite. For plausible values of c^, Vc\ is 
small enough even near the ISCO that the inspiral time of the 
clumps there is long enough to be consistent with the observed 
coherence of the kilohertz QPOs. 

3. RESULTS AND DISCUSSION 

The model developed here of the effect of the inward drift 
of clumps on the frequencies of the kilohertz QPOs has three 
parameters: the gravitational mass M and spin frequency I's 
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Fig. 2. — Comparison of fits of the SPBF gas dynamical model described in 
the text (solid curves) with the h.v-V2 correlations observed in four sources. 
The best-fitting parameter values ai'e Sco X-1: M = 1.59 A/q, = 352 Hz, 
cjv = 5 X 10-3, x^/dof = 95.5/46; 4U 1608-52: M = 1.76 Af©, 
= 343 Hz, Cjv = 2.8 x 10"^, x^/dof = 16.9/10; 4U 1728-34: 
M = 1.74:5 Mq, Cjv = 6 X 10""', x^/dof = 8.5/6; and 4U 1820-30: 
i^s = 279 Hz, CJV = 10^5, x^/dof = 7.0/18. For4U 1728-34, Us was set 
equal to the 364 Hz frequency of its burst oscillations; for 4U 1820—30, M 
was set equal to 2.0 Mq, the value indicated by the plateau in its frequency- 
flux relation (see text). 

of the neutron star and the coupling coefficient cjv ■ If M or i/g 
are known, there are only two parameters. Figure 2 shows the 
best fits of the model to the frequencies of the kilohertz QPOs 
in Sco X-1,4U 1608-52, 4U 1728-34, and 4U 1820-30. The 
model gives good fits (x^/dof of 0.4 to 2.1, not including sys- 
tematic errors), masses (1.59 to 2.0 A/q), spin rates (279 to 
364 Hz), and coupling coefficients for all four sources. 

The refined version of the SPBF model presented here is con- 
sistent with our previous prediction (MLP98) of a plateau in 
the kilohertz QPO frequency-X-ray luminosity relation when 
the sonic radius approaches the ISCO. In particular, the fit to 
the 4U 1820—30 data shown in Figure ^is consistent with the 
interpretation that the plateau in the QPO frequency vs. X-ray 
flux relation reported in this source (Zhang et al. 1998; Kaaret 
et al. 1999) indicates that (1) there is an ISCO around this neu- 



tron star and (2) the asymptotic frequency of the upper kilohertz 
QPO is the frequency of the ISCO. 

Finally, we note that the fundamental frequency in the SPBF 
model is the orbital frequency at the sonic radius, which is typ- 
ically ^ 1, 000 Hz. Hence, even the largest observed difference 
between and (70 Hz in 4U 1728-34) and the lai-gest 
observed variation of Ai^ (100 Hz in 4U 1608-52), although 
highly significant, represent only very small (5%-7%) devia- 
tions of the kilohertz QPO frequencies from the values com- 
puted in the circular orbit approximation used previously. 

In summary, we have shown that in the sonic -point beat- 
frequency model, the inward radial motion of the accretion 
flow causes the frequency separation Ai' between the kilohertz 
QPOs to be less than the spin frequency of the neutron star. 
Aiy is equal to I's only if the inflow velocity at the sonic radius 
is negligible. The inflow velocities expected at Rgp naturally 
produce the small (^0.2%-5%) corrections to the frequencies 
predicted for circular orbits (MLP98) needed to explain the dif- 
ferences between Ai/ and i's and the decreases of Ai/ with in- 
creasing 1^2 observed in some sources. We have developed a 
simple physical model of the effect of inflow on the frequencies 
of the upper and lower kilohertz QPOs. The model fits well the 
frequencies of the kilohertz QPOs in four Z and atoll sources 
where they have been fairly accurately measured and gives rea- 
sonable values for the masses and spin rates of the neutron stars 
in these sources. 

We emphasize that the frequencies of the kilohertz QPOs pre- 
dicted by the SPBF model depend sensitively on the proper- 
ties of the spacetime near the neutron star, including the grav- 
itomagnetic torque produced by the spinning neutron star and 
whether there is an innermost stable circular orbit. Hence, if 
the sonic-point model is shown to be correct, measurements 
of the frequencies of the kilohertz QPOs can be used not only 
to determine the properties of neutron stars but also to detect 
and measure gravitational effects in the strong-field regime (see 
MLP98). 
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